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Abstract 

We have obtained the exact asymptotics of the determinant 
deti<r,s<i [^("^^l^^) 6xp(i0) Jr,sj • Inverse symbolic computing methods 
were used to obtain exact analytical expressions for all terms up to rela- 
tive order L"^* to the leading term. This determinant is known to give 
weighted enumerations of cyclically symmetric plane partitions, weighted 
enumerations of certain families of vicious walkers and it has been conjec- 
tured to be proportional to the one point function of the 0(1) loop model 
on a cylinder of circumference L. We apply our result to the loop model 
and give exact expressions for the asymptotics of the average of the num- 
ber of loops surrounding a point and the fluctuation in this number. For 
the related bond percolation model at the critical point, we give exact ex- 
pressions for the asymptotics of the probability that a point is on a cluster 
that wraps around a cylinder of even circumference and the probability 
that a point is on a cluster spanning a cylinder of odd circumference. 

1 Introduction 

The binomial determinant 



det 

l<r,s<L 



r + s~2 
r- 1 



exp(i6') Sr 



(1) 



gives a weighted enumeration of certain types of nonintersecting lattice paths 
and cyclically symmetrical plane partitions 8J. Consider a family of noninter- 
secting lattice paths on an L x L square lattice subjected to the constraints: 

1. Paths are allowed to move one step to the right or downward. 

2. Paths are required to start from a point (0, k) with < fc < L — 1. 

3. A path starting at (0, k) is required to end at the point (fc, 0). 
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Figure 1: A typical configuration of nonintersecting lattice paths enumerated 
by the determinant in ([1} for L = 4. Paths starting at Ai have to end at Ei. 



See Fig.[l]for an example. The determinant ([T]) gives a weighted enumeration of 
all possible lattice paths satisfying the above constraints, where a configuration 
containing s paths is given a weight of exp[i0(L — s)]. 

Families of lattice paths of this type are in bijection with cyclically sym- 
metrical plane partitions [H]- A plane partition of an integer N is an ar- 
ray of nonnegative integers nj^k, such that rijj^i^k < ?^j.fe, nj,k+i < nj,k and 
^ — Sjli SfcLi '^i.fc- Plane partitions can be represented by a pile of unit 
cubes by introducing x,y,z coordinates in and placing at position (j, k, 0) 
a stack of rij^k unit cubes. A cyclically symmetric plane partition is a plane 
partition whose representation as a pile of cubes is symmetric under a cyclic 
permutation of the x,y,z coordinates. The bijection maps lattice paths of the 
above type to cyclically symmetric plane partitions that fit in an L x L x L box. 
Families of lattice paths containing n paths are mapped to cyclically symmetric 
plane partitions with n cubes on the main diagonal. 

The determinant ^ has been evaluated exactly for the cases a multiple 
of 7r/3 13 EH]- For general 9, an approximate asymptotic expansion was 
obtained by Mitra and Nienhuis [Hj. They numerically studied the dense 0(1) 
loop model [7] on a cylinder. The dense 0(1) loop model on square lattices can 
be defined as follows. At each vertex, the four edges meeting there are connected 
with equal probability in either of the two ways shown in Fig.O The set of edges 
connected to each other will then form closed loops with probability 1. In Fig. 





Figure 2: The two vertices of the dense 0(1) loop model 
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Opart of a typical example of a loop configuration is shown on a 6 x cx3 lattice. 
Here we have imposed periodic boundary conditions in the horizontal direction 
which gives the lattice the topology of a cylinder with a circumference of 6. If 
the circumference is chosen to be odd, then there will be "loop" spanning the 
entire length of the cylinder. 




Figure 3: Part of a typical configuration of the 0(1) loop model on a 6 x oo 
cylinder, periodic in the horizontal direction and extending to infinity in both 
vertical directions. 



Mitra and Nienhuis conjectured an exact formula for the probability P{L, m) 
that a point on a cylinder with circumference L is surrounded by m loops [18j . 
From that conjecture they were able to relate the generating function for the 
probability distribution to the determinant ([1]) [19]. Defining 



L/2 



(P{L, 9)=J2 ^(^' 2" cos™(6') 



and 



D{L,0) = exp{-ieL/2) det 

l<r,s<L 



r + s-2 
r- 1 



exp(i0) 5r 



the conjecture is 
Conjecture 1 



D{L, 



Aht{LY 



if L is even, 



D(L, 



2cos(f) AHTiLf 



if L is odd. 



(2) 
(3) 



(4) 
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Here A}n{L) is the number of L x L half-turn symmetric alternating sign ma- 
trices [111 mug: 

. ^ VT'3(fc-l)!fc!(3fe-l)!(3fc + 2)! ^ . 

Aht(L) = 2 ^ — — 5 — - if L IS even, 

4(2fc- l)!^(2fc + l)P 

^22 ""^^ 

Am{L) = ff ^^^1' if L is odd. 

f=\ my' 

We note that the function D{L,9) defined in Eq. [3] is an even function of 6. 
This follows from a symmetry property of the coefficients of the characteristic 
polynomial of the Pascal matrix. If we put 



det 

l<r,s<L 



r + s- 2 
r - I 



— xS. 



r.s 



Y.''kx\ (6) 



then = (—1)'^'^ C2,/2-fc, see [T7]. For reference later on in this paper, we 

note that the total number of n x n alternating sign matrices, A[n)^ is given by 

[II [12 

io (' + ")' 

The asymptotic behavior of correlations in the 0(1) loop model can be cal- 
culated using nonrigorous Coulomb gas techniques, see [IHllinilH]. As shown 
in [IH], this leads to the following conjecture. 

Conjecture 2 



n— — oo 



where the functions Cn{L, 6), which are not the same for odd and even L, have 
an asymptotic expansion in powers of l/L, such that the term proportional to 
l/L vanishes. 



Conjectures [Hand [2] imply for ~tt < 9 < tt: 

(9) 



D{L, 9) = A{9) ( ^ 1 L^/36-3«V(4.^) \i + o( _-J__^ 

^ ' ' '\ A ) [ \L'^inl2,3{l--\e\/7T)] J 

In Rcf. the amplitude A{9) was evaluated numerically. It was observed 
that: 

Conjecture 3 The amplitude A{9) in the expansion of D{L, 9) is the same 
function for odd and even L. 

Exact values for A{9) were obtained for the cases 9 a multiple of 7r/3 using the 
known exact evaluations of D{L,9). 

The results of this paper were obtained by making two additional assump- 
tions. 
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Assumption 1 The functions C'n{L,6) in ^ are analytic in 9. 

This was also mentioned in [19] . This assumption, combined with the periodicity 
of 0(L, 9) imphes that the C„(i, 9) are analytic continuations of each other 

C„(L,6') =Co(L, 61 + 27rn). (10) 

Another assumption we make is 

Assumption 2 The asymptotic expansion oj Cq{L,9) in powers of l/L con- 
tains only even powers of l/L. 

This is supported by the known exact evaluations of D{L, 9) for 9 a multiple of 
7r/3 and by numerical computations which we discuss in the next section. 

Using Conjectures [1] and [5] and Assumptions [T] and [21 we were able to gen- 
erate very accurate numerical data for the terms in the asymptotic expansion 
of Co{L,9). This allowed us to guess the exact expressions in the asymptotic 
expansion up to relative order L~^^ to the leading term. The asymptotic ex- 
pression for D{L,9) is given for both even and odd L by 

oo 

D{L,9)^ ^ (-l)"V(L,0 + 27rn), (11) 

n— — OO 

where 

/,M)^(|)*«P,-5.,-i,)(5^)^*-l*o(:.A).(,_^) 

(12) 

Here G{z) is the Barnes G-function and the R2k{S) are polynomials 

of degree 2fc -f 2. For fc = 1, . . . , 7 they are given by Eq. (jHJ in Section [31 For 
odd L this formula is implied by the validity of this formula for even L and 
Conjecture [31 

In the remainder of this paper we will discuss in detail how the above ex- 
pression was obtained. In Section [21 we explain how we obtained the numerical 
data for the terms in the asymptotic expansion. In SectionlSlwe discuss guessing 
exact expressions from the numerical data and how that led to our result. In 
Section [H we apply the obtained asymptotic expansion to the loop model and 
compute the asymptotics for a number of correlations in this model. 



2 Generating high precision numerical data 

For convenience, we will first only consider the case of even L. We will return 
to the case of odd L later on in this paper. From ([3]), (O, and (|T]) it follows that 

^MlL3eV(4.^)-i/i2 ^ g f;Afe(^? + 27rn)i-3"'-3"'^/-^ (13) 

HTV / n=-oo k=0 
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Table 1: Values for D{9) for 9 a multiple of 7r/3. The expressfons for A{L) and 
Awt[L) are given in d?]) and ^ 
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even L 


odd L 





Aht(2L) 


Aht(2-L) 


7r/3 


AHT(i)' 




27r/3 




ML) 


TT 








where the Ak{9) are even functions of 9. To compute these functions, we gen- 
erated the characteristic polynomial of the Pascal matrix up to L = 200, which 
took a few days on an ordinary PC. Using these polynomials we could quickly 
evaluate the left-hand side of the above expression. This allowed us to extract 
the Ak{9) and their first few derivatives. 

Let us first see what we can conclude from the exact expressions for D{L, 9 — n^) , 
see Table [H We can expand these expressions in powers of 1/L by rewriting the 
product of factorials in terms of Barnes G- functions [H [5J Ul [S] [5] and using the 
asymptotic expansion of the Barnes G-function, see e.g. Ref. [T^ for an efficient 
method to do this. We note that the asymptotic expansion of the Barnes G- 
function given in some papers contains a sign error, in particular in Refs. Plll9j. 
The correct expansion is 

log[G(z + 1)] Q log(^) - i) + ^ log(2^) ^ - ^ log(^) + C'(-l) 

R (14) 

\ ^ i>2fc+2 

^ 4fc(fc 1) z2fc ■ 

Here the Bk are the Bernoulli numbers, and C'(^) is the derivative of the zeta 
function. This can be derived using the methods presented in [2| . A simple sign 
error leading to an incorrect minus sign in front of the summation over k was 
made in that article. 

It turns out that the asymptotic expansions of the expressions in Table [T] 
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contain only even powers of \/L: 
31/12 



0(L,O) 



■exp(-C'(-l)) 



25/36 

13294838545991999 



1 



127 



2041055 



5184L2 53747712L4 
645434518069131955571 



107538127903 
835884417024L6 



17332899271409664L8 89853749822987698176Lio 
272944577297197688875376083 26385460676926169502575757887765 



2794811034494209364066304L12 

TT 

3 

2tt 



14488300402817981343319719936L1 



= 1 



L 



231/ 



36 



TT 



exp(-C'(-l)) 



3 J ' 35/12 p(^/3)2 
225307455655 1215802858094435 



7 



23983 



16317695 



576L2 663552X4 127401984L6 



293534171136L8 169075682574336Lio 
19038476800109154745 204450994938396835527815 



194775186325635072L1 



28/3 



TT 



3 r(i/3r 

45379702784 



112190507323565801472L14 
8 464 228352 



77553152 



81L2 6561L4 1594323L6 129140163L8 
122234658136064 156017791843041280 



10460353203L 
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2541865828329X12 2058911320946492.14 



(15) 



The absence of odd powers of 1/L suggests that for general 9, Ak{9) = for 
odd k. Our numerical results are consistent with this (the fact that AiiO) — 
is a Coulomb gas prediction). From (fTS]) we see that the A2k{d + 2-Kn) will also 
give rise to odd powers of 1/L when 9 — Q or 6 = 27r/3 for odd n: 



A2fc[27r(2n+1)] =0, 



A2k 



27r 

— + 27r(2n+l) 



0. 



(16) 



From the leading terms of the above expansions we can read-off ylo(^p) for 
p = 0, 1, 2, 3. Note that Ao(7r) is half the leading term of (t){L, tt) because n = 
and n = —1 contribute equally to the summation in (jl3p when 9 = n. We can 
read-off some more Ak (f p) for fc > 0, but it is not possible to find more special 
cases of the form Ao(^p -t- 2?™) due to mixing with the Ak{^p). To calculate 
these values, we must lift the degeneracy between these terms, which can be 
done by considering the derivatives of (|13p w.r.t. 9. 
Our strategy to find Ao{9) was as follows: 

obtain accurate fits for special values of 



1. Using (derivatives of) Eq. ((T 
Ao{9) and its derivatives. 

2. Try to guess exact expressions for these approximate values. 

3. Try to get information about the zeroes of Ao{9) from the fits. 

4. Assuming that Aq{9) is of finite order, any guess about the location of all 
the zeroes (and their multiplicities) corresponds, by Hadamard's factor- 
ization theorem [T^, to a guess for ^40(6*) up to a factor exp(polynomial). 
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Such a guess can be easily checked by evaluating summations over the ze- 
roes of the form J2n "^n/ (o^n ~ where the q;„ are the zeroes of Ao{9), 
the m„ the multiplicities of the zeroes, and p is an integer larger than 
the order of Ao{9). Such summations can be expressed in terms of the 
derivatives of Ao{9) at 6* = a as follows. The contour integral 

^o(^) dz ^^^^ 



c{B) Ao{z){z-af' 

where C{R) is a circle with the origin as its center and radius R, will tend 
to zero for R — > oo if p is chosen larger than the order of A{6). It then 
follows from the residue theorem that 



m„ 1 dp 



(a„ - af (p- 1)! dzP 



log[Ao(z)] 



(18) 



If, using cither numerical or exact expressions for the derivatives of Aq[9)^ 
we verify the validity of the above formula for some particular a and p, 
we can construct a candidate function for ^o(^) smd verify that it fits the 
data. Note that it is theoretically possible for the last test to fail. E.g., it 
could be that p is not larger than the order of ^o(^) and the guess about 
the zeroes is also wrong in such a way that (fT5|) is satisfied by accident. 

5. If (HI fails, we could try to apply some transformation to Aq{9), e.g., by 
subtracting some function from it and consider the zeroes of the trans- 
formed function. Or we could abandon the attempt to find the function 
from its zeroes altogether and do some sort of brute force search instead. 



3 Finding exact expressions from numerical data 

An important result we obtained was the derivative of 740(6*) at = vr/S. The 
numerical value we obtained by fitting the derivative of ()13|) was 



A'^^i^ « -0.31250232645180558377609866618987578223847934983479 . . . 

(19) 

We conservatively estimated the accuracy of this number to be ^ 30 digits. To 
get an idea what the analytical expression of this number could possibly be like, 
let us consider the derivative of the function <j>{L,9). It follows from fS]) that 
the leading asymptotics of the derivative is 

*'('"l)--^»(l)-2;'°«<^)' (»' 

Now, for finite L the value 0'(i, ^) is just an algebraic number. The fact that 
on the right-hand side there is a logarithmic term in L strongly suggests that 
the "constant term" ^o(-|) contains a term proportional to Euler's constant 
with the same amplitude as the coefficient in front of the log(L). If we subtract 
this contribution of Euler's constant from our fit to ^o(f ) '^^ S^^ 

^of^l+TT- ~ -0.22063560015265159339645643211799769082485006927415 . . . 

(21) 
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It is not so clear what this number is, but perhaps we should multiply it by tt 
because the term log(L) + 7 appears with a factor 1/{2tt): 



27r. 

This looks familiar: 



- 0.69314718055994530941723212145817656757963215693787 . 



-log(2) = - 0.69314718055994530941723212145817656807550013436026. 

(22) 

So, we arrive at the following conjecture: 



For general 9: 



ct^'{L,e)^A',{e)-^A,{e)\og{L). (24) 



And we thus expect that A'f^{9) / Aq{9) contains a term — 30/(277^) 7. After 
subtracting this contribution we may not be able to recognize the remaining 
term as easily as in the case of 6* ~ 7r/3. The problem we need to solve is in 
general: We suspect that some number x which we know to limited precision 
is a rational linear combination of some constants yi . . . We want to know 
if this is true, in which case we want to find the linear combination. This type 
of problem can be solved using so-called integer relation algorithms. We used 
the Lenstra-Lenstra-Lovasz lattice reduction algorithm (LLL algorithm) [TB] for 
this purpose, implemented in Mathematica by the function LatticeReduce. In 
Appendix [X] we give a brief description of this algorithm and explain how it can 
be used as an integer relation algorithm. 

Using the LLL algorithm were able to find more exact values for the loga- 
rithmic derivative of Ao(0), such as 

^o(f ) _ -27 + 31og(3) 

A'^{^) ^ ~37-61og(2) + 31og(3) 
Ao(^) 27r 

The pattern of zeroes was not very clear at this stage. Fitting to the second 
derivative of (fT3|l dX 9 — Q yielded very small coefficients for the log(i)^ 
terms, consistent with these terms being exactly zero. This suggests that there 
are zeroes at 27rn for n 0. Fitting to the first and second derivative of (fT5)) 
at 6* = 27r/3 yielded very small coefficients for the log(L) L"*' and log(L)^ L^^ 
terms, respectively. This then suggests additional zeroes at 27rri ± 27r/3. Note 
that for odd n this coincides with 

As we can see from ([T51) . the Ak{9) rapidly diverge as a function of k. This 
makes it difficult to extract accurate values for Ao{9) for 9 > 3tt. We could only 
verify that the derivatives at 2tt, 47r/3 and 87r/3 were nonzero. Clearly it is not 
possible for all the zeroes to be simple zeroes because then Ao{9) would have 
to be a trigonometric function up to a polynomial factor, which is inconsistent 
with the results we found so far. E.g. Aq{9) cannot be a product of functions like 
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r(a + b6) T{a — h6). A product of functions of the form T{a + bO'^) is consistent 
with the derivatives we found, but is not consistent with the equal spacing of 
the zeroes. A product of Barnes G-functions is a more plausible possibility, 
however, then the multiplicities of the zeroes must increase linearly, but that is 
not what we see at the first few zeroes. 

We proceeded with attempting to find expressions for higher derivatives. The 
presence of Euler's constant in the first logarithmic derivative suggests that in 
higher (logarithmic) derivatives we can expect to find polygamma functions. 
Using the LLL algorithm we found 



Am 

Ao{0) 



3 37 31og(3) , ^1(3) ^i(i) 



27r2 27r2 



27r2 



+ 



67r2 



67r2 



(26) 



where ipp{z) is the polygamma function of order p. If the order of Ao{9) is 2 or 
higher then we need to consider higher derivatives to compute summations over 
the zeroes. If the order of Ao{9) is less than 4, then ^„ mn/a-n is given by —4 
times the coefiicient of 9'^ in the series expansion of log[Ao(6')] around ^ = 0. 
After some trial and error we found using the LLL algorithm the expression: 



1 ^ 
4! ^6*4 



log[Ao(0)] 



1 



6>=0 



576774 



V'3 



V'3 



27 
32^ 



C(3) . (27) 



We can reproduce this result if we assume that all the zeroes of {9) are at 27rn 
and 2im± 27r/3 and have a multiplicity of |n|. The first three positive zeroes at 
47r/3, 27r and 87r/3 then have a multiplicity of 1 as we observed earlier. To see 
that this works, let us write the summation over the zeroes as 



E 



m„ _ 1 
1 



1 



Ei X ^ n -i ^ n 



.n=l 



00 ^ 00 ^ 00 ^ 

~^ / 1 \ ^ 5^ / . 1 



247r4 



„=i(n-i) n=o{n+\) _ 



y 1 -y^— 

.„=! \n - 3) „=o + 3) 



(28) 



We can simplify this expression further using 

E ^ + E '( ITS + E '( ~7\3 = E 7;;T3 = 27C(3) , 
and the identity: 



1 

1 

in A- 



which allows us to write 

00 



„=l("-3) „=o("+3) 



6"H3 



and 



(29) 



(30) 



(31) 



(32) 



10 



These simplifications yield the expression: 

1 

STr^''"' ' 1447r4 



(33) 



We see from ([27|) that this is indeed the same as —4 times the coefficient of 9'^ 

of \og[Ao{e)]. 

The function G(l + z) has zeroes at the negative integers z = —n with 
multiplicity n and is of second order. The conjecture about the zeroes of ^0(6*) 
and the assumption that the order of Aq{6) is less than four is thus equivalent 
to the statement: There exist constants a and h such that 



Ao(0)=exp(a + 6e^)G(l + -)G(l--)G(- + - 



^G[- + — ]G[-~ — 



27r 



27r 



27r 



27r 



We verified that this fits the data exactly. We found 



G -- — 
3 27r 



(34) 



Ao{e) =2"/363i/36^i/3 



exp 



?C'(-i))r(^ 



27r 



G 1 + — G — 



27r 



27r 



27r 



xG^ + — G--:^ G- + :^ G--:^ 



27r 



27r 



(35) 



Next, we turned our attention to the A2k{S)- We found that the ratios 
A2k{6) IAq{6) appeared to be polynomials of degree 4fc. We observed that the 
summation over k in (|13p can be simplified as follows. If we define the polyno- 
mials Q2k{d) as 



Ao{9)exp 



.k=Q 



12k 



{0) 



I2k 



E 

fe=0 



A2k{d) 



(36) 



then Q2k{d) is of degree 2fc-|-2. Our results so far can be summarized as follows. 
For even L, the function 0(L, 9) is given by an asymptotic expansion of the form: 



E 



g2fc(6' + 27rn) 



L2fe 



(37) 



What about odd LI In Ref. [19] it was observed that the coefficient of the 

^_3eV(47r^) + l/12 ^^^^ ^ Ag<id(0) = 2c^e/2) M0)- An 

equation of the same form as ((57)) for odd L would thus imply that there exist 
polynomials Q§^'^(^) such that 

%/3 



</.(L,0) 



2cos( 



^ (-1)" Ao(i, + 27rn) i-3(e+2n.)V(4.^)+i/i2 



n— — 00 



X exp 



E' 



27m) 



(38) 



for odd i. From (jj]) we see that for even L, the function D{L,9) is obtained 
by multiplying (f>(L,9) by ylHT(i) , while for odd L we must multiply by 
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2cos(|) v4ht(-^)^- As was observed in [12], this leads to the same resuh for 
the leading term of D{L,9) because A^n^L) for odd L is 3"^/'* times AhtC^) 
for even L in the large L limit. Our results therefore imply that the asymptotics 
for D{L,0) for both even and odd L are given by the formula: 

oo 

^ (-l)"V(i,e + 27rn) (39) 



where 
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f(L. ,) ^ (I) " e.p(-5<'(-l)) j + 2^ j - 2, 



(40) 

The polynomials R2k{0) for = 1, . . . , 7 are given by 



R2{0) - 
ReiO) = 
Rio{0) = 

Rl2{0) = 

i?i4(e) 



77 76*2 lir 



15552 1447r2 647r4 

245 157612 296*4 1816*6 



559872 129607r2 11527r4 12807r6 
1103 1349612 3599614 ggg^e 32756*^ 



40310784 2449447r2 311047r4 69127r6 143367r8 

793135 116807612 1010096*4 474796*^ 431716i« 616216*1" 



4353564672 20995207r2 2799367r4 6220807r6 368647r8 1228807rio 
93651593 37400096*2 18680836*^ 3010916*^5 18585136** 



130606940160 103926247r2 13996807r4 9331277^ 2764807r8 
12397736*1° 11841716*12 



1843207rio 9011207ri2 
2884889645 680610916016*2 1100185696*^ 7548141430^ 



940369969152 232133428807r2 226748167r4 167961607r6 
4548716216*8 9316522936*1" 311937316*12 230575816*1^ 



104509447r8 99532807ri" 8847367ri2 59637767ri4 
2213492219141 24715026056*2 830194155316*4 308696349196**^ 
135413275557888 ~ 765275047r2 ~ Tl42567040^ 45349632 tt^ 
101009167736*8 969362374916*1° 356196713896*12 1052933156*14 



447897607r8 627056647rio 398131207ri2 5898247ri4 

4530052916*16 



367001607ri6 



(41) 



We note that the summation in the exponent in (|40|) is a divergent asymp- 
totic expansion of some unknown function. If wc truncate this expansion at 
some finite order, then the summation over n in (|39p will diverge. However, 
since the summation over n converges very fast (independently of i), this is not 
an issue when doing numerical computations. 
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4 Application to the 0(1) loop model and bond 
percolation 

In this section, we consider the leading asymptotics of some correlations in the 
0(1) loop model and bond percolation model on the cylinder that can now be 
evaluated exactly. The bond percolation model at the critical point is related to 
the dense 0(1) loop model via the bijection shown in Fig. H) We can see from 




Figure 4: The mapping of a loop configuration to a bond configuration of the 
corresponding bond percolation problem. A bond is either put on an edge of 
the square lattice formed by the • or on the dual edge orthogonal to it on the 
square lattice formed by the o. 

this figure that the loops of the dense 0(1) loop model correspond to bound- 
aries of percolation clusters. From ^ we can see that the generating function 
(/)(L, 9 = 7r/2) equals the probability that a point in the loop model is not sur- 
rounded by any loops. This probability thus corresponds to the probability in 
the percolation model that a point is on a cluster that wraps around the cylin- 
der for even L, while for odd L it gives the probability that the point is on a 
cluster that spans the cylinder. This probability was given in [18] for even L as 
P « 0.81099753L-5/48[i + 0{L-^/^)]. Using §7^) we can evaluate the leading 
asymptotics as 



35/48 y/fXT/T) 



(42) 



It follows from ((57|) that the leading asymptotics for the probability that a point 
is on a cluster that spans the cylinder is •\/3/2 times the above expression. 

Next, we consider the number of loops surrounding a point. By repeatedly 
differentiating both sides of Eq. ^ at 9 = 7r/3, we can express the moments 
of P{L,m) in terms of derivatives of (j){L,9) a,t 9 — tt/3. These derivatives can 
be evaluated using the derivatives of ((57)) and (|38p at 6* = tt/3. Denoting the 
number of loops surrounding a point by N, we have 

{N) = -i^[7 + log(4L)] [1 + 0{L-')] for even L, 
(TV) = |-- + + log(4L)]| [1 + 0{L-')] for odd L. 

Using the first and second derivatives at 9 — tt/3 the fluctuation in N can be 
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calculated. This is given by 

l + log(3) , 7Ai(l/6) 



{N')-iNf={-l- 



27r2 187r2 
[1 + 0{L-^)] for even L, 

l + log(3) , Vi(l/6) 



27r2 187r2 
[1 + 0(1.-2)] for odd L. 



3V37r 27r2 



[7 + log(4L)] 



[7 + log(4L)] 



(44) 



5 Discussion 

We have obtained the exact asymptotics of the characteristic polynomial of 
the Pascal matrix. The form of the asymptotics was partially known due to 
a conjecture relating D{L,9) to (j){L,6). This was used to obtain accurate 
numerical data for the amplitude of the leading power of L. The LLL lattice 
reduction algorithm was used to guess exact expressions for special values of 
the amplitude which ultimately enabled us to find exact expressions for the 
amplitude for all 9 and the amplitudes of the subleading terms up to relative 
order L~^^ to the leading term. 

We have applied this result to the 0(1) loop model and the related bond 
percolation model on the cylinder. Asymptotic expressions for the average of 
the number of loops surrounding a point and the fluctuation in this number 
were obtained. We also obtained the leading asymptotics of the probability 
that a point is not surrounded by any loops, which for even L corresponds to 
the probability that a point is on a cluster that wraps around the cylinder, 
while for odd L it gives the probability that a point is on a cluster that spans 
the length of the cylinder. 

The methods presented here can likely be used to obtain exact asymptotic 
expressions for a large class of binomial determinants. E.g., in [5], exact expres- 



for arbitrary in- 



sions for the determinant deti<r,s<L (''^^,^1 ^) + exp(i6') (5^,, 
teger m and 9 a multiple of 7r/3 were obtained. The asymptotics for general 9 
appears to have a similar structure as for the case m = investigated in this 
article. 
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A Solving integer relation problems using the 
LLL algorithm 

The LLL algorithm solves the following problem. Suppose we are given n linear 
independent vectors vi . . . v„ in R*^ with d > n. Consider the set of all integer 
linear combinations of the v, . Such a set is called a lattice and the v,- are a 
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basis of the lattice. The basis of a lattice is not unique in general. Given the 
Vj", the LLL algorithm finds short, approximately orthogonal, basis vectors for 
the lattice spanned by the Vj . 

In an integer relation problem, one is given a real number x which one wishes 
to express in terms of real constants yk '■ 



E 



rkVk, (45) 



where the are rational numbers. In our case, x is known to limited precision 
while the yk are known to arbitrary precision. We want to know if for given x 
and yk such a relation exists, in which case we want to find the r^. The LLL 
algorithm can be used to solve this problem as follows. We first multiply our 
number x and the constants yi ■ ■ .yn by some power of ten and take the floor 
to get integers: 

x' = Lio%j . ^ ' 

We choose p such that x' contains only those digits of x that we know to be 
correct. We then define the basis vectors: 

= + y ■e„+2 for I < i < n, 
Vn+i = e„+i + X e„+2- 

where is the rth unit vector of R"+^. We take these n+1 vectors vi . . . v„_|_i 
as the input of the LLL algorithm. The output will be a set of short basis 
vectors for the lattice spanned by the v^. If there is a short basis vector b in 
the output for which the {n + 2)nd component is 0(1) in n and the (n + l)st 
component is nonzero, then we have very likely found an exact relation between 
the yi and x: 

n 

^bkyk + bn+ix = 0. (48) 

k=l 

The larger the number n of constants is, the larger we need to choose p to be 
able to detect a relation. 

The LLL algorithm is not the best available integer relation algorithm. The 
PSLQ algorithm [TT] is faster and will yield lower bounds on the sizes of the 
coefficients if no relation is found. However, unlike the PSLQ algorithm, the 
LLL algorithm can also be used to find an unknown function given the action of 
a finite number of linear functionals. This would be useful in the event we need 
to move on to point [5] of the program described on p. [71 Suppose an unknown 
function / : 5*1 — > 5*2 for arbitrary vector spaces Si and 5*2 satisfies the relations: 

Lrf = Xr (49) 

for I < r < m where the L^. are linear functionals and the Xr are real numbers. 
If / is suspected to be a rational linear combination of functions gk'. 

n 

f^^rkgk, (50) 

k=l 
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where n > m, we can use the LLL algorithm to find the if they exist. We 
define 

yk,r = Lrgk, 

= [WPyk,r\ , (51) 

for 1 < fc < n and 1 < r < m, where p is limited by the accuracy of the Xr- We 
then take as input for the LLL algorithm the basis vectors defined as 

m 

v; = e., + ^ J/- .^e„+i+,. for 1 < z < n, 

(52) 

If a short output basis vector b has a nonzero {n + l)st component and the 
components (n + 1 + r) for r — 1, . . . , m are all 0(1) in n, then that vector very 
likely defines an exact relation between / and the gk- 

n 

Y,bk9k + bn+if^0, (53) 

k=l 

provided p was chosen large enough. 
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